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ABSTRACT 

We present the results of a ground-based, high spatial resolution infrared 18 /jm imaging study of 
nearby luminous infrared galaxies (LIRGs), using the Subaru 8.2-m and Gemini South 8.1-m tele- 
scopes. The diffraction-limited images routinely achieved with these telescopes in the Q-band (17-23 
/xm) allow us to investigate the detailed spatial distribution of infrared emission in these LIRGs. We 
then investigate whether the emission surface brightnesses are modest, as observed in starbursts, or are 
so high that luminous active galactic nuclei (AGNs; high emission surface brightness energy sources) 
are indicated. The sample consists of 18 luminous buried AGN candidates and starburst-classified 
LIRGs identified in earlier infrared spectroscopy. We find that the infrared 18 /im emission from the 
buried AGN candidates is generally compact, and the estimated emission surface brightnesses are 
high, sometimes exceeding the maximum value observed in and theoretically predicted for a starburst 
phenomenon. The starburst-classified LIRGs usually display spatially extended 18 /im emission and 
the estimated emission surface brightnesses are modest, within the range sustained by a starburst 
phenomenon. The general agreement between infrared spectroscopic and imaging energy diagnostic 
methods suggests that both are useful tools for understanding the hidden energy sources of the dusty 
LIRG population. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — galaxies: starburst - 
infrared: galaxies 



1. INTRODUCTION 

A large number of bright galaxies that emit most 
of their radiation in the infrared (8-1000 /im; Ltr, > 
IO^Lq), peaking at ~60 /im, were found with the IRAS 
infrared sky survey (Soifer et al. 1987). These are called 
luminous infrared galaxies (LIRGs) or ultraluminous in- 
frared galaxies (ULIRGs) if the infrared luminosity ex- 
ceeds Lir > 10 12 L Q (Sanders et al. 1988a; Sanders & 
Mirabel 1996). The strong infrared emission means that 
powerful energy sources are present, hidden behind dust, 
which absorbs most of the primary energetic radiation; 
the heated dust grains re-emit this energy as infrared 
thermal radiation. These energy sources can be star- 
bursts (energy release by nuclear fusion reactions in- 
side stars), or activity within an active galactic nucleus 
(AGN; energy conversion from gravitational energy gen- 
erated by a mass-accreting supermassive black hole to 
radiative energy), or some combination of the two. Once 
the primary energetic radiation is converted to heat, 
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identifying the original energy source becomes an obser- 
vational challenge. 

An effective method of investigating the energy sources 
of LIRGs is to measure their emission surface bright- 
nesses. In a starburst, the radiative energy generation 
efficiency of nuclear fusion is only ~0.5% of Mc 2 , where 
M is the mass of material used in the nuclear fusion re- 
action, and c is the speed of light. The emission surface 
brightness of a starburst is modest and cannot exceed a 
certain threshold (~10 13 Lq kpc -2 ), as determined obser- 
vationally (Werner et al. 1976; Meurer et al. 1997; Soifer 
et al. 2000; Riechers et al. 2009) and predicted theoret- 
ically (Elmegreen 1999; Thompson et al. 2005; Younger 
et al. 2008). In contrast, the radiative energy generation 
efficiency of an AGN is much higher, at 6-42% of Mc 2 , 
where M is the mass of accreting material (Bardeen 1970; 
Thorne 1974). In an AGN, high luminosity can be gener- 
ated from a very compact area around a mass-accreting 
supermassive black hole, producing a very high emission 
surface brightness with >10 13 L e kpc -2 . 

For the dusty LIRG population, high spatial resolu- 
tion imaging observations in the infrared are a suitable 
choice to constrain the emission surface brightnesses of 
energy sources, because dust absorbs most of the primary 
energetic radiation, infrared dust re-radiation dominates 
the observed spectral energy distribution, and so traces 
the intrinsic energetic radiation luminosities in LIRGs. 
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For a galaxy with given infrared luminosity, when the 
emission is spatially unresolved, a higher spatial resolu- 
tion image can constrain the size of the emission region 
more strongly, and can provide a more stringent lower 
limit for the emission surface brightness of an energy 
source. Ground-based 8-m class telescopes are better 
suited for this purpose than space-based infrared satel- 
lites with small apertures, as the former provides bet- 
ter spatial resolution than the latter. To investigate the 
spatial distribution of infrared emission of LIRGs from 
the ground, A-band (8-13 /xm) and Q-band (17-23 /xm) 
Earth atmospheric windows are possible choices. Higher 
sensitivity can be achieved in the A-band than Q-band, 
due to a smaller background level from Earth's atmo- 
sphere. Thus, A-band observations are more suitable 
for probing the overall morphology of infrared emission 
of LIRGs, including both nuclear compact and diffuse 
extended components, to a fainter flux level (Soifer et 
al. 2001; Alonso-Hcrrero et al. 2006; Diaz-Santos et al. 
2008; Siebenmorgen et al. 2008). However, for the pur- 
pose of constraining the emission surface brightnesses 
of energy sources in LIRG's nuclei, we believe that Q- 
band observations are more powerful, for the following 
reasons. First, the long wavelengths of the Q-band are 
significantly less affected by turbulence in the Earth's at- 
mosphere than shorter- wavelength bands. Thus, in the 
Q-band, diffraction-limited image sizes (~0'.'5) are read- 
ily achievable using ground-based 8-m class telescopes 
under even modestly good weather conditions. Hence, 
the point spread function (PSF) at Q is stable over time, 
and the investigation of the intrinsic spatial extent of 
compact emission in LIRGs is very reliable. Second, in 
the iV-band, (1) the contribution from poly cyclic aro- 
matic hydrocarbons (PAH) emission features, originat- 
ing in star-forming activity, can be strong, relative to 
dust continuum emission, and (2) infrared emission pow- 
ered by an obscured energy source can be highly flux- 
attenuated by the strong 9.7 /xm silicate dust absorption 
feature (Armus et al. 2007; Imanishi et al. 2007a; Iman- 
ishi 2009; Veilleux et al. 2009; Imanishi et al. 2010a). 
Suppose the case that ^90% of the total infrared (8-1000 
/xm) luminosity of a LIRG is dominated by a compact, 
highly-obscured buried AGN, with a small (^10%) con- 
tribution from spatially-extended star-forming activity. 
While the A-band emission from the buried AGN can 
be highly flux-attenuated, spatially-extended PAH emis- 
sion is unattenuated. If we observe such a LIRG in the 
A-band, the fraction of the spatially-extended emission 
component, relative to the compact one, can be substan- 
tially higher than the actual energetic contribution to the 
total infrared luminosity, possibly misidentifying such a 
source as starburst-important. This possible ambiguity 
is largely overcome, if we observe the LIRG in the Q- 
band, because (1) the contribution from PAH emission 
features is much reduced, compared to the A-band, and 
(2) the 18 /xm silicate dust absorption feature is signifi- 
cantly weaker than the 9.7 /xm feature. We thus choose 
the Q-band for our infrared imaging observations. 

With the advent of high-sensitivity Q-band imaging 
cameras attached to 8-m class telescopes at good astro- 
nomical observing sites, Q-band high spatial resolution 
imaging observations of LIRGs have now become tech- 
nically feasible. We have therefore performed such ob- 
servations, using the Subaru 8.2-m telescope located on 



Mauna Kea, Hawaii, and the Gemini South 8.1-m tele- 
scope located at Cerro Pachon, Chile. Throughout this 
paper, H = 75 km s _1 Mpc -1 , = 0.3, and 17a = 0.7 
are adopted. 

2. TARGETS 

Due to the strong Earth atmospheric background emis- 
sion in the Q-band (17-23 /xm), sensitivity is still lim- 
ited even with 8-m class telescopes. Therefore, we chose 
nearby LIRGs that are bright in the Q-band as our main 
targets. As LIRGs have highly concentrated nuclear gas 
and dust (Sanders & Mirabel 1996; Soifer et al. 2000), 
many of the putative AGNs are likely to be obscured 
along virtually all directions at the inner regions and 
deeply buried. Such buried AGNs lack the narrow line re- 
gions at the 10-1000 pc scale, photoionized by the AGN 
radiation, and so optical spectroscopy fails to discover 
AGN (Seyfert) signatures. Our primary scientific goal is 
to search for such optically elusive buried AGNs in op- 
tically non-Seyfert LIRGs, through the detection of high 
emission surface brightness energy sources. In fact, high 
spatial resolution Q-band imaging, as a tool for quan- 
titatively estimating the emission surface brightnesses 
of energy sources, is particularly effective at picking up 
buried AGNs with a dust covering factor of about unity, 
because almost all primary energetic radiation is con- 
verted to dust re-radiation, which can be probed with 
infrared observations. We selected LIRGs with luminous 
buried AGN signatures in X-ray, infrared, or millimeter 
spectra, as well as starburst-classified LIRGs (without 
buried AGN signatures) and optical Seyfert 2s (= ob- 
scured by torus-shaped dusty medium) of interest for 
comparison. A few very nearby well-studied galaxies 
with Lir < 10 11 Lq were also included in our sample, be- 
cause their proximity allows us to investigate the infrared 
emission morphology in more detail on a physical scale. 
Table 1 tabulates the 18 observed objects. The sample is 
heterogeneous. Our scientific goals are (1) to strengthen 
any buried AGN signatures suggested in previously ob- 
tained data, and (2) to determine if the emission sur- 
face brightnesses of energy sources in starburst-classified 
galaxies are indeed below the maximum limit sustained 
by a starburst phenomenon (^1O 13 L kpc -2 ), as widely 
argued. 

3. OBSERVATIONS AND DATA ANALYSIS 

3.1. Q-band Observations 

Gemini Q-band (17-23 /xm) imaging observations were 
made using T-ReCS (Thermal-Region Camera Spectro- 
graph; Telesco et al. 1998) on the Gemini South 8.1- 
m telescope, at Cerro Pachon, Chile. T-ReCS uses a 
Raytheon 320 x 240 pixel Si:As IBC array, with a pixel 
scale of 0'.'09. The field-of-view is 28'.' 8 x21'.'6. The detec- 
tor was read out in correlated quadruple sampling (CQS) 
mode (Sako et al. 2003). The Qa filter (18.3 /xm; 17.6- 
19.1 /tm) was used to trace the ^18 /xm emission from 
LIRGs. The standard chop-nod technique (telescope 
nodding and secondary mirror chopping) was employed 
to subtract the background emission from Earth's atmo- 
sphere and the telescope. Object signals were taken at 
four different positions: (Nod, Chop) = (1, A), (1, B), (2, 
A), and (2, B). However, for T-ReCS data, only one chop 
position was guided, and compact, diffraction-limited- 
images were obtained. At another chop position, the 
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telescope was unguidcd, and stellar images were clearly 
degraded and spatially extended. Only the data taken 
at the guided chop position were used for our analysis, 
to investigate the intrinsic spatial distribution of the in- 
frared 18 /im emission in a reliable manner. The chop 
and nod throw was 15". For some sources with spatially 
extended structures, the position angle of T-ReCS was 
arranged to simultaneously cover the emission from these 
interesting regions. The frame rate was 25.8 ms, and to- 
tal on-source exposure time was 15-31 minutes. Table 2 
summarizes our detailed observing log. Weather condi- 
tions were very good (low background level) in 2008 and 
2009, but during the 2010 observing runs, a high Earth 
atmospheric background level allowed us to observe only 
bright sources. 

Subaru Q-band imaging observations were made using 
COMICS (Kataza et al. 2000; Okamoto et al. 2003) on 
the Subaru 8.2-m telescope (lye et al. 2004) atop Mauna 
Kea, Hawaii, under clear weather conditions. The details 
are shown in Table 2. The precipitable water vapor value 
was as low as <1 mm during the 2008 observing run, 
and 3-4 mm during the 2009 run. The Q17.7 (17.7 ^m; 
17.25-18.15 fim) filter was used to study the ^18 ^m 
emission from LIRGs. COMICS employs a 320 x 240 
Si:As IBC detector for imaging observations. The pixel 
scale is 0'.'13, providing a field of view of 42" x 32" , if the 
full array is used. However, as the background emission 
was large for Q17.7-band imaging observations, we used 
a partial-readout mode, to avoid saturating the array. 
Only 50-120 rows were read out, depending on weather 
conditions, particularly on the precipitable water vapor 
value. The standard chop-nod technique was employed 
to subtract background emission. The frame rate was 16- 
31 ms, depending on the actual background level at the 
time of observations, and the total on-source exposure 
time ranged from 16 to 60 minutes. 

Our Gemini T-RcCS data were reduced using the 
IRAF-based analysis tool provided by the Gemini obser- 
vatory 4 . After inspecting individual frames to confirm 
that no bad frames were included, the task "mireduce" 
was used to combine frames to improve the signal-to- 
noise ratios. For some frames, a stripe-pattern noise was 
recognizable. We used the task "miclean" to correct for 
this noise pattern. 

The Subaru COMICS data were analyzed using the 
IRAF-based software package provided by the COMICS 
team 5 . First, a flat- field image was created from frames 
with no object signals, and then individual frames were 
divided by the flat image. Bad pixels and pixels hit by 
cosmic rays were identified, and signal values in these pix- 
els were interpolated from the signals of the surrounding 
pixels. The frames were summed to improve the signal- 
to-noise ratios. 

For both Gemini South T-ReCS and Subaru COMICS 
data, flux calibrations were made by comparing the sig- 
nals of the targets with those of nearby bright standard 
stars compiled by Cohen 6 . The estimated flux levels 
of the target objects with sufficient signal at different 
nod/chop positions agreed to within <0.1 mag. Thus, 

4 http://www.gemini.edu/sciops/instruments/t-recs/data- 
format-and-rcduction?q=nodc/10145 

5 http:/ /canadia. ir.isas.ac.jp/comics/open/rbin/rbin. html 

6 http://www.gcmini.edu/sciops/instrumcnts/mir/Cohen_list.html 



the adopted choice for the absolute flux level of stan- 
dard stars may represent the largest uncertainty when 
compared with other flux measurements reported in the 
literature. 

The observed nearby bright standard stars were also 
used to estimate the point spread function (PSF) at the 
observed wavelengths, at the time close to the target ob- 
servations, in a sky direction similar to the targets. The 
standard star images in the Qa-band (18.3 /jm) taken 
with Gemini South T-ReCS generally showed diffrac- 
tion ring patterns around bright central peaks. Those 
in the Q17.7-band taken with Subaru COMICS often 
showed the so-called Mitsubishi-pattern (three diamonds 
around central bright cores) . Both of these observed fea- 
tures strongly suggest that nearly diffraction-limited im- 
age sizes were routinely achieved at ~18 /im for both the 
Gemini South and Subaru telescopes. 

Standard stars were observed before and after the tar- 
gets. The FWHM values differed slightly between the 
two datasets. Although the effects of seeing (Earth at- 
mospheric turbulence) at ^18 /im are much smaller than 
at shorter wavelengths, all images can obviously suffer 
from a non-negligible amount of seeing-originated PSF 
broadening. Thus, it is reasonable to assume that the 
smallest FWHM is the most representative of the intrin- 
sic PSF (Radomski et al. 2008). As standard stars are 
bright, we can estimate the FWHM by dividing the com- 
bined final frame into individual frames. We adopted the 
smallest FWHM values found in some individual frames 
as the PSF at ~18 //m (FWHM PSF ). 

The PSFs of standard stars were estimated using the 
IRAF task imexamine, a. Following Radomski et al. 
(2008), the Moffat function (Moffat 1969) was basically 
adopted, but direct FWHM estimates were sometimes 
considered as a consistency check. In general, both val- 
ues agreed for the data with sufficiently strong signals. 

For Gemini South T-ReCS data, Radomski et al. 
(2008) measured the FWHM = 0'.'53 ± 0'.'039 at the same 
Qa-band. Our adopted FWHM values for the standard 
star (FWHMpsf) ranged from 0'.'515 to 0'.'565, consis- 
tent with the measurements by Radomski et al. (2008). 
The Superantcnnae galaxy is at declination —72°, and is 
very far south even for the Gemini South telescope site, 
resulting in a larger air mass, even at the meridian, than 
other LIRGs. The FWHM of HR 7383, the standard 
star of Superantennae, was actually slightly larger than 
the FWHMs of other standard stars, possibly because of 
the telescope primary mirror misalignment caused by a 
larger air mass. For Subaru COMICS data, the adopted 
FWHM values for the standard star PSFs fell in the 
0'.'41-0"45 range. Thus, it is important to use the FWHM 
of a standard star observed at an air mass similar to the 
air mass of each target. 

As our target objects, LIRGs, are generally fainter 
than bright standard stars, it was difficult to measure the 
FWHM values in individual frames in a reliable manner, 
due to the low signal-to-noise ratios. Thus, the FWHM 
values were reliably measurable only in the combined fi- 
nal frames. For the Gemini South T-ReCS LIRG data, 
one final dataset consisted of a ~15-min combined expo- 
sure, and one or two such datasets were taken. For mod- 
erately bright LIRGs, we inspected individual frames in 
each dataset, but no obvious continuous image drift was 
detected. Therefore, we measured the FWHM values in 
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the final dataset. For objects with two such datasets, a 
smaller FWHM value was adopted. 

For Subaru COMICS data of LIRGs, one dataset con- 
sisted of 4-minute exposures, in which we could identify 
object signals and roughly measure the object peak po- 
sitions. For some LIRGs, the peak position continuously 
shifted along one direction at a maximum distance of a 
few pixels. Since the air mass changed as observations 
proceeded, this could be explained by the differences in 
refraction between the optical (guiding camera) and mid- 
infrared light (COMICS camera) due to the Earth's at- 
mosphere. We corrected for this shift manually at both 
the pixel scale (using the IRAF task imshift) and the 
sub-pixel scale (using the IRAF task imshift and ima- 
lign). We then produced a final combined frame, and 
measured the FWHM values. The sub-pixel scale shift 
often provided larger FWHM values than the pixel scale 
shift, possibly because of the signal split into surrounding 
pixels during the sub-pixelization process. We adopted 
the smaller FWHM values as the ~18 /im emission sizes 
of LIRGs (FWHM L i R g). 

The intrinsic 18 fim emission sizes of LIRGs 
(FWHMjntrmsic) were estimated using the formula 
FWHM intrinsic 2 = FWHM LmG 2 - FWHM PSF 2 . While 
the FWHM values for the standard stars were derived 
from individual short exposure frames, those for LIRGs 
were obtained from combined frames with longer expo- 
sure times. The FWHM values in the final combined 
frames are inevitably larger than the smallest FWHM 
values found in some individual short exposure frames, 
due to possible telescope tracking error (drift) and/or 
possible inclusion of poorer seeing data in some frac- 
tion of individual frames than the remaining majority of 
frames. Thus, the FWHMi ntr insic 2 values for LIRGs es- 
timated in this way are conservative upper limits, or the 
estimated emission surface brightnesses (= flux divided 
by size) of LIRGs are conservative lower limits. How- 
ever, our main conclusions are robust to this ambiguity 
because the presence of a luminous AGN is considered 
only if the lower limit of the emission surface brightness 
of the energy source in the LIRG substantially exceeds 
the maximum allowed value for a starburst phenomenon 
(see §5). 

Finally, as mentioned by Radomski et al. (2008), spec- 
tral energy distributions differ between LIRGs and stan- 
dard stars, in such a way that standard stars are gener- 
ally bluer than LIRGs. Hence, the effective wavelength in 
the observed filters can be slightly different (i.e., shorter 
for standard stars), possibly providing different PSFs be- 
tween LIRGs and standard stars. However, the widths 
of the Qa (Gemini South T-ReCS) and Q17.7 (Subaru 
COMICS) filters used are so narrow that this effect is 
practically negligible (Radomski et al. 2008). Even if it 
were non- negligible, this effect would result in a slightly 
larger PSF for LIRGs (red) than bluer standard stars, 
providing upper limits (lower limits) for the emission size 
(emission surface brightnesses) of LIRGs. Again, this 
ambiguity does not alter our conclusions (see §5). 

3.2. 10 fim Observations 

When precipitable water vapor values were high (>3 
mm) during Subaru COMICS observing runs, we per- 
formed TV-band (8-13 (im) imaging observations of bright 
sources, mostly unobscured quasars (QSOs), as a backup 



program. The A8.8 (8.4-9.2 fim) or Nil. 7 (11.2-12.1 
(im) filters were used, because meaningful data could be 
obtained for bright sources in these filters (compared to 
the Q17.7 filter) even under high precipitable water con- 
ditions. It is interesting to measure the emission surface 
brightnesses of these QSOs (i.e., known luminous unob- 
scured AGNs), based on the ~10 fim observations. Table 
3 shows our observing log. 

For the 2006 observing run, only 160 rows were read 
out. The standard chop- nod technique was employed to 
subtract the background emission. The exposure times 
were 41-100 ms, and the total on-source exposure time 
ranged from 13 to 53 minutes. For the 2009 observing 
run of the LIRG Mrk 266, observations were made only 
at the nod-1 position, using the full array, to simultane- 
ously cover both the northern and southern nuclei with 
a separation of ~10". 

The FWHMs of standard stars and target objects were 
derived in the same way as outlined for the 18 /im obser- 
vations. The IRAF task imexamine was used, and the 
Moffat function was basically assumed. For bright stan- 
dard stars, data were divided into individual frames, and 
the smallest FWHM value was adopted as the intrinsic 
PSF because any unexpected factors (tracking errors or 
sudden seeing degradation during parts of the exposure 
time) work to increase the FWHM values. As the QSOs 
were relatively bright, FWHM values were derivable in 
individual frames (100-s exposure). The smallest FWHM 
value was adopted. As the exposure time of individual 
frames was longer for QSOs (100 s) than standard stars 
(5 s), this procedure still provides conservative upper lim- 
its for the intrinsic spatial extent of ~10 fim emission of 
QSOs. The ~10 fim emission of QSOs was dominated 
by compact components. The LIRG Mrk 266 displayed 
a spatially extended emission, compared to the adopted 
PSF. 

3.3. Spitzer IRS Spectroscopy 

The presence of spatially extended, faint fuzzy emis- 
sion, undetectable with our ground-based 18 /im imag- 
ing observations, can be investigated by comparing our 
data with the spectra taken by Spitzer IRS LL2 (14.0- 
21.3 fim) (Houck et al. 2004), if available. The aperture 
size of IRS LL2 with 10'.'5 is sufficiently large to cover 
any extended diffuse infrared 18 fim emission. As our 
targets are well-studied LIRGs, Spitzer IRS LL2 spectra 
at 18 fim are available for many of them (Brandl et al. 
2006; Roussel et al. 2006; Armus et al. 2007; Imanishi 
et al. 2007a; Wu et al. 2009). However, to our knowl- 
edge, the flux-calibrated Spitzer IRS LL2 spectra have 
not been published for IRAS 20551-4250 and VV 114 
E, although LL2 observations have been performed ac- 
cording to the Spitzer observing log. We have analyzed 
archival spectra of these sources. For IRAS 15250+3609 
and NGC 1377, strong 18 /im silicate dust absorption 
features were found in the flux-calibrated LL2 spectra 
(Armus et al. 2007; Roussel et al. 2006). We re-analyzed 
the IRS spectra to estimate the optical depth of the 18 
fim silicate dust absorption feature based on our method 
(Imanishi et al. 2007a; Imanishi 2009; Imanishi et al. 
2010a), which is required for dust extinction correction 
to the AGN emission (§5.1). For NGC 5253, Spitzer 
IRS SH high-resolution (R ~ 600) spectrum at 18 fim, 
with spectral mapping mode (total covered area is 22'.'2 
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x 33'.'4) is available (Beirao et al. 2006). This is used 
for the flux comparison at 18 /an. For ESO 602— G025, 
NGC 7592, and Arp 193, wc find no record of Spitzer 
IRS spectroscopic observations at ~18 /zm. 

The Spitzer IRS observing log of the sources analyzed 
in this study is shown in Table 4. All four modules, SL2 
(5.2-7.7 Atm), SL1 (7.4-14.5 /zm), LL2 (14.0-21.3 /zm), 
and LL1 (19.5-38.0 /zm) were used. We analyzed all 
data to obtain the full 5-35 /zm low-resolution (R ~ 100) 
spectra. 

Data analysis was performed in a standard manner, 
using methods similar to those employed by Imanishi et 
al. (2007a), Imanishi (2009), and Imanishi et al. (2010a). 
The latest pipeline-processed data products available at 
the time of our analysis were used. Frames taken at po- 
sition A were subtracted from those taken at position 
B to remove background emission, consisting mostly of 
zodiacal light. The spectra extracted for the A and B 
positions were then summed. Wavelength and flux cali- 
brations were made on the basis of the Spitzer pipeline- 
processed data. For the SL1 spectra, data at A Q b s > 14.5 
/zm in the observed frame are invalid (Infrared Spectro- 
graph Data Handbook, version 1.0), and were therefore 
removed. For the LL1 spectra, data at A bs > 35 /zm 
were not used due to excessive noise levels. These data 
are not necessary for our scientific discussion. When flux 
discrepancies between SL1 and LL2 were discernible in 
such a way that the LL2 flux (1075 wide slit) was larger 
than the SL1 flux (377), we adjusted the SL1 (and SL2) 
flux to match the LL2 flux. 

4. RESULTS 

4.1. Imaging 

Figures 1 and 2 show the infrared 18 /zm images of 
LIRGs obtained with Gemini-South T-RcCS Qa-band 
(18.3 /zm) and Subaru COMICS Q17.7-band (17.7 /zm) 
observations, respectively. Photometric measurements 
are summarized in Table 5. 

Most of the bright standard stars show clear signatures 
of a diffraction ring or partial ring, strongly suggesting 
that diffraction-limited images were routinely obtained 
at 18 /zm under the normally good, frequently achiev- 
able weather conditions at the Gemini-South and Subaru 
telescopes. Hence, at 18 /zm, the PSF is stable, and the 
possible spatial extents of LIRGs can be investigated in a 
reliable manner, with the data only minimally affected by 
the time variation of the PSF. For LIRGs, only VV 114 
E, ESO 602-G025, NGC 4945, NGC 1614, and Arp 193 
display clearly detectable spatially extended structures. 
The 18 /zm emission of NGC 5253 is spatially compact 
but resolvable. The remaining sources are dominated by 
spatially compact emission, which is not clearly distin- 
guishable from the PSFs of standard stars. 

For VV 114 E, the two faint knots below the E-2 nu- 
cleus in Figure 1 were not recognizable in the low spa- 
tial resolution 15 /zm image taken with ISO, but similar 
structures were seen at 1.1 /zm (Le Floc'h et al. 2002), 
suggesting that they are starburst knots situated between 
the VV 114 E and W merging nuclei. 

For NGC 4945, the elongated structures along the NE 
to SW direction were also seen at shorter infrared wave- 
lengths (Moorwood et al. 1996; Marconi et al. 2000). 

4.2. Spectroscopy 



Figure 3 shows the final Spitzer IRS spectra. IRAS 
15250+3609, IRAS 20511-4250, and NGC 1377 show 
clearly detectable 18 /zm silicate dust absorption features, 
while the feature is weak in VV 114 E. The estimated op- 
tical depths, based on our method (Imanishi et al. 2007a) 
are n 8 - 0.6, 0.7, and 0.85 for IRAS 15250+3609, IRAS 
20511-4250, and NGC 1377, respectively. 

5. DISCUSSION 

5.1. Intrinsic Emission Surface Brightnesses of Energy 

Sources 

5.1.1. Comparison of the 18 /zm Flux Levels 

Our ground-based 18 /zm imaging observations may 
not be particularly sensitive to spatially extended dif- 
fuse emission outside the detected compact cores. To 
test the presence of a possible fuzzy component, we com- 
pare our 18 /zm fluxes for the core components with 18 
/zm fluxes measured with Spitzer IRS spectra (1075 aper- 
ture or even larger for spectral mapping mode; §3.3) in 
Table 5. Given the ~20% flux uncertainties of Spitzer 
spectra (Spitzer Infrared Spectrograph Data Handbook, 
version 1.0), wc see significantly (>>20%) smaller fluxes 
in our measurements only in VV 114 E and NGC 4945, 
both of which clearly display spatially extended emission 
structures. Therefore, we find no strong evidence that 
the contribution from spatially extended diffuse emission 
outside the compact core component is important, but is 
simply missed with our ground-based imaging observa- 
tions, in the majority of the observed LIRGs. 

For IRAS 08572+3915 and 15250+3609, our fluxes are 
slightly (~100 mJy) smaller than Spitzer IRS measure- 
ments, although the discrepancy is close to the possible 
flux calibration uncertainty of Spitzer IRS. Even if this 
residual flux with ~100 mJy at 18 /zm were due to spa- 
tially extended diffuse starburst emission below our de- 
tection limit, the expected 60 /zm flux is ~1 Jy, if wc 
adopt typical spectral energy distributions of starburst 
galaxies at 17-30 /zm (Brandl et al. 2006), and IRAS 
60-/zm to 25-/zm flux ratios for the prototypical star- 
burst galaxy, M82 (~4). The expected 60 /zm flux for 
the possible extended starburst component is a factor 
of >7 smaller than the measured IRAS 60 /zm fluxes of 
IRAS 08572+3915 and 15250+3609 (Table 1), suggesting 
that any undetected spatially extended starburst emis- 
sion cannot dominate the large infrared luminosities of 
these LIRGs. 

For VV 114 E, the sum of the 18 /zm fluxes for 
the compact components (E-l and E-2) measured with 
our ground-based observations are >30% lower than the 
Spitzer IRS measurement. Thus, a substantial amount 
of diffuse, faint emission component may be present in 
VV 114 E. 

For NGC 4945, the nuclear 18 /zm flux (178 x 178 
region of the brightest point A in Figure 1) is ^200 mJy. 
When signals in the region of discernible emission (879 x 
376) are integrated, the 18 /zm flux becomes ~1000 mJy. 
These fluxes are still smaller than the Spitzer IRS 18 /zm 
flux of ^1500 mJy (1075 aperture), suggesting that low 
surface brightness diffuse emission may exist, and which 
is not detectable with our ground-based 18 /zm imaging 
observations in NGC 4945. 

5.1.2. Size and Emission Surface Brightness 
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The sizes of the 18 /im emission in apparent and phys- 
ical scales, as estimated by the method described in §3.1, 
are summarized in Table 6. For LIRGs whose core emis- 
sion is clearly extended (VV 114 E, ESO 602-G025, 
NGC 5253, NGC 1614, and Arp 193), the actual val- 
ues are adopted for the emission size. However, for the 
remaining LIRGs, image sizes are only marginally larger 
than the PSFs of the corresponding standard stars. Part 
of this effect may be artificially induced by the longer ex- 
posures time of LIRGs compared to standard stars (§3.1). 
For LIRGs whose 18 /im image sizes are not clearly dis- 
tinguishable from the PSFs, upper limits are shown for 
the sizes of the 18 /im emission regions. The emission 
surface brightnesses are estimated by dividing the 18 /im 
continuum luminosity (y\j v ) with the physical scales of 
the 18 /im emitting regions (Table 6). 

For sources for which only A-band imaging observa- 
tions were performed under high prccipitable water va- 
por conditions (§3.2), we estimated the sizes of the 10 
/im emission regions and emission surface brightnesses 
from the A-band data, using methods similar to those 
outlined above. Table 7 summarizes the results. 

5.1.3. Flux- attenuation Correction of 18 [im Continuum 

Emission 

For LIRGs with clear 18 /im silicate dust absorption 
features, it is suggested that the 17-18 /im continuum 
emission in the rest-frame (outside the 18 /im silicate 
dust absorption feature) is obscured and flux-attenuated 
by foreground absorbing dust. Proper correction of this 
flux attenuation is necessary to estimate the intrinsic, 
extinction-corrected luminosities of the dust-obscured 
energy sources from the 18 /im observations. 

As a dust extinction curve as long as 18 /im is not well- 
constrained observationally, we combine the theoretical 
dust extinction curve at A > 7 /im with the observed 
extinction curve at shorter wavelengths. We adopt the 
following relations. 

1. Tis/Tg.7 ~ 0.3 (Chiar & Tielens 2006; Imanishi et 
al. 2007a) 

2. t 9 7 /A v = 0.06 (Roche & Aitken 1984, 1985; Rieke 
& Lebofsky 1985) 

3. A 7Aim /A v ~ 0.0225 (Draine & Lee 1984; Lutz et al. 
1996; Nishiyama et al. 2008, 2009) 

4. A 17fim /A 7lJ , m - 1.03 (Draine & Lee 1984) 

From these relations, we obtain Ai7 Mm = 1.3 x ris. 
For example, for a source with ti 8 = 1, the 17-18 /im 
continuum flux is attenuated by a factor of ^3.3. By 
adopting the T\s values based on our estimate (Imanishi 
et al. 2007a; § 4.2 of this paper), the corrected emission 
surface brightnesses are derivable for LIRGs that exhibit 
clear 18 /im silicate dust absorption features. The cor- 
rected values are shown in Table 6. 

5.1.4. Dust Covering Factor Correction 

Infrared observations can probe only emission ab- 
sorbed by dust. For buried AGNs embedded in dust 
in virtually all directions, or starbursts whose individ- 
ual stellar energy sources are surrounded by local dust 



with a covering factor close to unity, the intrinsic lumi- 
nosities of primary energetic radiation can reasonably be 
estimated from the infrared 18 /im observations. How- 
ever, for AGNs with a dust covering factor significantly 
smaller than unity, a large amount of energetic radiation 
escapes without being absorbed by dust. Such unab- 
sorbed emission cannot be probed through infrared ob- 
servations. PG QSOs are unobscured Seyfert 1 AGNs, 
meaning that at least the direction along our line-of-sight 
is dust-free, in which case infrared observations can pro- 
vide only lower limits for the intrinsic emission surface 
brightnesses of energy sources. Following Risaliti et al. 
(2010), we apply a factor of 5 correction for unobscured 
Seyfert 1 AGNs. Namely, the intrinsic emission surface 
brightnesses of the energy sources are 5 times larger than 
those estimated from infrared 10 /im observations. These 
corrected values are also summarized in Table 7. 

5.1.5. Do the 18 /im Observations Precisely Trace the 
Spatial Distribution of 60 \im Emission 1 ? 

LIRGs usually show an infrared emission peak at ~60 
/tm, and their total infrared luminosities are dominated 
by emission from big dust grains in thermal equilibrium 
with radiation (Sanders & Mirabel 1996). In the Galac- 
tic diffuse inter-stellar medium, emission from PAHs and 
very small grains, which are transiently-heated and not 
in thermal equilibrium, contributes importantly to the 
infrared 18 /zm radiation (Desert et al. 1990). Then, a 
question may arise whether our 18 /im observations in- 
deed probe the spatial distribution of dust in thermal 
equilibrium heated by the dominant energy sources in 
LIRGs. We, however, believe that the 18 /im observa- 
tions are a good method to investigate the spatial dis- 
tribution of 60 /im emission in LIRG's nuclei, for the 
following reasons. 

First, in the case of a normal starburst, it is very 
likely that PAHs, very small grains, and big dust grains 
are spatially well mixed with stellar energy sources, and 
are heated and/or excited by the stellar radiation field. 
The spatial distribution of the emission from PAHs, very 
small grains, and big dust grains at LIRG's cores should 
not change dramatically on a few 100 pc scale. 

Next, in the case of an AGN, PAHs are largely de- 
stroyed by strong X-rays, so that PAH contribution to 
the observed infrared radiation becomes negligible (Voit 
1992). Furthermore, high radiation density around an 
AGN (1) heats big dust grains in thermal equilibrium to 
high temperature, and (2) destroys very small grains se- 
lectively, both of which make the contribution from big 
dust grains dominant to the infrared 18 (im radiation. 

Finally, in a buried AGN, the energy source (= a mass- 
accreting supermassive blackhole) is more centrally- 
concentrated than the surrounding dust, so that dust 
temperature is higher at the inner part, closer to the 
central energy source, than the outer part. Assuming ap- 
proximately blackbody radiation, emission at a shorter 
infrared wavelength comes from more inside than that at 
a longer infrared wavelength (Figure 2 of Imanishi et al. 
2007a). Namely, the 18 /im emission size is intrinsically 
smaller than the 60 fim emission size, around a buried 
AGN with centrally-concentrated energy source geome- 
try. The same situation could occur, if an exceptionally 
centrally-concentrated, extreme starburst is considered 
(Figure le of Imanishi et al. 2007a). However, even for 
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a luminous energy source with 10 L Q , the intrinsic size 
difference between the 18 /zm and 60 /zm emitting regions 
should appear on a physical scale much smaller than 100 
pc, and should not be discernible in our spatial resolu- 
tion of a few 100 pc physical scale for the majority of the 
observed LIRGs. 

Summarizing, the spatial extent of the 18 /zm emission 
is taken to be representative of that of the 60 /zm emission 
at LIRG's cores on our physical resolution scale. The 
emission surface brightnesses estimated from our 18 /zm 
imaging observations can put important constraints on 
the nature of energy sources at LIRG's cores. 

5.2. Individual Sources 

In Tables 6 and 7, sources are categorized into three 
groups: (1) The lower limit for the emission surface 
brightness of the energy source is >10 13 L Q kpc -2 , the 
maximum value achieved by a starburst phenomenon. 
(2) The actual value is substantially smaller than the 
starburst upper limit. (3) The lower limit is smaller than 
the starburst upper limit. We discuss these results on an 
individual source basis, by comparison with other data 
reported in the literature. 

5.2.1. Nearby Sources with Energetically Important 

AGNs 

The following five sources at z < 0.1 are diag- 
nosed to be AGN-important through infrared 2.5-35 
/zm spectroscopy: NGC 1377, IRAS 08572+3915, IRAS 
15250+3609, Superantennae, and IRAS 20551-4250 
(Table 1, column 10). The estimated emission surface 
brightnesses of the energy sources exceed ~10 13 L© kpc -2 
(the maximum value of a starburst phenomenon) by a 
factor of >3 for IRAS 08572+3915 and 15250+3609 (Ta- 
ble 6), supporting the presence of luminous buried AGNs 
in these sources. For Superantennae, IRAS 20551—4250, 
and NGC 1377, the lower limits of the emission surface 
brightnesses only barely exceed ~10 13 L Q kpc -2 . As Su- 
perantennae is classified optically as a Seyfert 2 (Table 
1), the dust covering factor around that AGN could be 
smaller than the other buried AGN candidates (IRAS 
20551-4250 and NGC 1377), so that the intrinsic emis- 
sion surface brightness could be higher than our estimate 
in Table 6. In these three sources, although the presence 
of a luminous AGN is a possibility, the identity of the 
energy sources is not strongly constrained from our 18 
/zm imaging observations. 

5.2.2. Distant Sources with Energetically Important 

AGNs' 

The five ULIRGs, IRAS 00091-0736, 00188-0856, 
01004-2237, 01298-0744, and 04103-2838 are also clas- 
sified into buried-AGN-important galaxies through in- 
frared 2.5-35 /zm spectroscopy, but their redshifts are 
larger than z > 0.1 (Table 1). The lower limits of the 
emission surface brightnesses are below or only barely 
larger than ^10 13 L Q kpc -2 , so that the constraints on 
the energy sources based on our infrared imaging meth- 
ods are weak. Observational constraints on the size of 
emission regions are practically limited by the apparent 
spatial extent of the infrared emission, which corresponds 
to a larger physical scale for more distant sources. Higher 
spatial resolution imaging observations at 18 /zm, using 



a proposed 30-m class extremely large telescope (Toku- 
naga et al. 2010; Okamoto et al. 2010) will be required 
to better constrain the energy sources of ULIRGs at z > 

0. 1. 

5.2.3. Unobscured AGNs 

In Table 7, the emission surface brightnesses of five 
PG QSOs (= unobscured Seyfert 1 AGNs) are estimated 
from 10 /zm imaging data, and their lower limits are close 
to or barely exceed ^1O 13 L kpc -2 . Although these re- 
sults are consistent with the AGN-dominatcd picture of 
these PG QSOs, further quantitative discussions are dif- 
ficult because of (1) the uncertainty of the dust covering 
factor correction, and (2) the possibly lower PSF stabil- 
ity of the 10 /zm imaging data compared to our 18 /zm 
results (§1). 

5.2.4. Starburst and AGN Co-existing Sources 

In four LIRGs, NGC 2623, NGC 7592 W, VV 114 E, 
and Mrk 266 N, the presence of an obscured AGN is 
suggested, but starburst activity is also strong (Table 1, 
column 10). The estimated lower limits for the emis- 
sion surface brightnesses are far below ~10 13 L Q kpc -2 , 
and so our data do not provide evidence for energetically 
important AGNs. 

NGC 4945 also belongs to this class. Due to its prox- 
imity (distance ~ 4.0 Mpc), the nuclear region is spa- 
tially resolved at small physical scales (1" corresponds 
to ^20 pc). The bright region, denoted "A" in Figure 

1, is observed as an elongated peak in the infrared L- 
band (3.5 /zm) and A-band (10 /zm) (Moorwood et al. 
1996; Marconi et al. 2000), but at 18 /zm, it is barely re- 
solved into two distinct peaks with a separation of ~0 / /5, 
possibly because of reduced dust extinction at 18 /zm. 
The left (NE) peak is even brighter; this might be the 
location of the putative buried AGN suggested from X- 
ray observations at >10 keV (Iwasawa et al. 1993; Done 
et al. 1996; Guainazzi et al. 2000). Two fainter distinct 
emission components are seen at the SW of nucleus A. 
A similar pattern is seen in the high spatial resolution 
infrared AT-band (2.2 /zm) image (Marconi et al. 2000). 
This is likely due to starburst knots. 

5.2.5. Starburst- dominant Sources 

Three LIRGs, NGC 1614, Arp 193, and ESO 
602— G025, have no obvious luminous AGN signatures 
in any observational data obtained to date, and are 
consistent with a starburst-dominant picture. These 
LIRGs were observed as a control sample to determine 
whether these starburst-classified galaxies show spatially 
extended structures and emission surface brightnesses 
within the starburst range. As expected, spatially ex- 
tended emission structures and low emission surface 
brightnesses are confirmed. Thus, for these LIRGs, our 
18 /zm imaging data support a starburst-dominant na- 
ture. 

NGC 5253 is a very nearby galaxy (distance ~ 3.9 
Mpc), dominated by a nuclear super star cluster (Gor- 
jian et al. 2001; Vanzi & Sauvage 2004; Alonso-Hcrrero et 
al. 2004). We observed this galaxy to determine whether 
emission surface brightness can substantially exceed the 
maximum value for a starburst phenomenon (~1O 13 L0 
kpc -2 ), if only a super star cluster region is extracted. 
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The 18 fxm emission is spatially compact, but is resolved. 
The emission surface brightness estimated based on our 
18 /im image is 1.5 x 10 13 L Q kpc~ 2 , roughly compara- 
ble to the maximum allowed value by a starburst phe- 
nomenon. 

5.3. Energy Sources 

A natural explanation for the energy sources with 
emission surface brightnesses considerably higher than 
is an AGN, because (1) AGN activity 
can have high radiative energy generation efficiency (§1) 
and (2) even a super star cluster produces only ^10 13 Lq 
kpc- 2 (§5.2.5). In fact, high emission surface bright- 
nesses are preferentially indicated for LIRGs with lumi- 
nous AGN signatures in other observational data, such 
as infrared 2.5-35 /im spectroscopy and/or X-ray obser- 
vations at E > 2 keV (Table 1, column 10). 

Thompson (2009) argued that the emission sur- 
face brightness of a starburst phenomenon can exceed 
~10 13 L Q kpc -2 , if the average dust temperature is T > 
200K. Theoretically, the maximum surface brightness of 
a starburst phenomenon is determined by the balance 
of stellar radiation pressure to dust and gravitational 
boundaries (Thompson et al. 2005). Under high stel- 
lar radiation density, dust and gas (collisionally coupled 
with dust) are expelled away, and further star formation 
activity is suppressed (Thompson et al. 2005). At a high 
average dust temperature, T > 200K, some fraction of 
small, fragile dust can be sublimated, decreasing the to- 
tal dust opacity to stellar radiation, and possibly increas- 
ing the maximum allowed emission surface brightness of 
a starburst. If this were the case in many LIRGs, the 
constraint of emission surface brightnesses of >>10 13 L Q 
kpc~ 2 would not serve as a strong AGN signature. 

Observationally, the spectral energy distributions of lo- 
cal LIRGs usually show peaks at ~60 fim, indicating 
typical dust temperatures of T ~ 50 K (Clements et al. 
2010). For dust temperatures of T > 200 K, the infrared 
emission should peak at ~15 /im, but no such sources 
are found in our sample (Table 1). Riechers et al. (2009) 
also calculated that the maximum dust temperature for 
a starburst is T ~ 88K. At present, the proposed star- 
burst with T > 200 K is not strongly supported from the 
observational point of view. 

The only way to reconcile the T > 200 K starburst 
scenario and the observed infrared spectral energy dis- 
tribution in LIRGs is that the T > 200 K starbursts are 
occurring only in the centrally concentrated regions of 
LIRG's nuclei, and are surrounded by large amounts of 
nuclear dust (similar to buried AGN geometry; Figure le 
of Imanishi et al. 2007a). In this geometry, the original T 
> 200 K starburst emission is once absorbed by the sur- 
rounding nuclear dust, and the re-emitted infrared dust 
emission indicates cooler temperatures, as observed in ac- 
tual LIRGs. However, LIRG's nuclear regions are known 
to contain large amounts of high-density molecular gas, 
sufficient to trigger star formation (Gao & Solomon 2004; 
Imanishi et al. 2007b; Wilson et al. 2008; Imanishi et al. 
2009). Thus, it seems reasonable for star formation to 
occur in the large volumes of LIRG's nuclear regions oc- 
cupied by high-density gas, and it is not clear why stars 
are formed only in the centrally concentrated regions of 
LIRG's nuclei. Therefore, we consider that high emis- 
sion surface brightness energy sources at LIRG's cores 



are good signatures of AGNs, if not definitive AGN evi- 
dence. 

6. SUMMARY 

The results of ground-based, high spatial resolution 
infrared 18 /im imaging observations of nearby LIRGs 
were presented. The intrinsic size of the infrared 18 
/zm emitting regions was estimated by comparison with 
the PSFs of standard stars, using (close to) diffraction- 
limited images. Emission surface brightnesses were then 
constrained. We drew the following conclusions from our 
observations: 

1. LIRGs with luminous AGN signatures seen in pre- 
vious infrared spectroscopy and/or X-ray observa- 
tions at >2 keV tended to show compact 18 ^m 
emission, while starburst-classified LIRGs with no 
obvious AGN signatures in previous observations 
displayed spatially extended morphologies with low 
emission surface brightnesses. For some fraction of 
the former LIRGs at z < 0.1, the lower limits of the 
emission surface brightnesses were constrained to 
>>10 13 L Q kpc -2 , higher than the maximum value 
sustained by a starburst phenomenon, as argued 
theoretically and actually observed. A luminous 
AGN is a natural energy source that can produce 
such a high emission surface brightness. 

2. For distant LIRGs at z >0.1, classified as AGN 
important through previous observations, 18 ^m 
emission was generally spatially compact. How- 
ever, the lower limits of the emission surface bright- 
ness were below or only barely exceed the upper 
limit expected from a starburst phenomenon. We 
could not place strong constraints on the presence 
of luminous AGNs for these distant sources, simply 
because the spatial resolution of ground-based 8-m 
class telescopes is still insufficient at 18 fim. 

3. We estimated the emission surface brightness of the 
super star cluster in the very nearby galaxy NGC 
5253 to be ~10 13 L© kpc -2 . This did not support 
the possibility that the emission surface brightness 
could be >>10 13 L Q kpc" 2 if only a super star clus- 
ter region is extracted. 
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TABLE 1 

Observed LIRGs and their IRAS-baseb infrared emission properties 



Object 




Redsbift 


fl2 




f60 


fioo 


log L IR 


f25/f60 


Optical 


Buried 








(Jy) 


(Jy) 


(Jy) 


(Jy) 


L© 




Class 


AGN 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


IRAS 00091-0738 




0.118 


<0.07 


0.22 


2.63 


2.52 


12.2 


0.08 (C) 


HII 


a,b,c 


IRAS 00188-0856 




0.128 


<0.12 


0.37 


2.59 


3.40 


12.3 


0.14 (C) 


LI 


a,b,c,d 


IRAS 01004-2237 




0.118 


0.11 


0.66 


2.29 


1.79 


12.3 


0.29 (W) 


HII 


a,b,c 


IRAS 01298-0744 




0.136 


<0.12 


0.19 


2.47 


2.08 


12.3 


0.08 (C) 


HII 


a,b,c 


IRAS 04103-2838 




0.118 


0.08 


0.54 


1.82 


1.71 


12.2 


0.30 (W) 


LI 


a,b,e 


irtAo Uoo/z+oylo 




U.Uoo 


U.oZ 


1 7(1 


*7 A Q 


4.oy 


1 O 1 


U.zo j 


T T 


a,b,c,d,f,g,h,i,j 


IRAS 15250+3609 




0.055 


0.16 


1.31 


7.10 


5.93 


12.0 


0.18 (C) 


LI 


b,c,h,k 


Superantennae (IRAS 19254- 


7924) 


0.062 


0.22 


1.24 


5.48 


5.79 


12.1 


0.23 (W) 


Sy2 


j,l,m,n 


IRAS 20551-4250 




0.043 


0.28 


1.91 


12.78 


9.95 


12.0 


0.15 (C) 


LI (HII) 


c,k,o,p,q 


NGC 1377 




0.006 


0.56 


1.93 


7.43 


5.95 


10.1 


0.26 (W) 


Unc 


r,s,t 


NGC 1614 




0.016 


1.38 


7.50 


32.12 


34.32 


11.6 


0.23 (W) 


HII 




NGC 2623 




0.018 


0.21 


1.81 


23.74 


25.88 


11.5 


0.08 (C) 


Unc 


u,v 


NGC 7592 




0.024 


0.26 


0.97 


8.05 


10.58 


11.3 


0.12 (C) 


Sy2 (W),HII (E) 


k 


Arp 193 




0.023 


0.25 


1.42 


17.04 


24.38 


11.6 


0.08 (C) 


LI 




ESO 602-G025 (IRAS 22287 


-1917) 


0.025 


0.27 


0.91 


5.42 


9.64 


11.3 


0.17 (C) 
0.16 (C) A 


LI 




VV 114 E 


0.020 


1.03 A 


3.65 A 


22.93 A 


31.55 A 


11.7 A 


HII 


i,w 


NGC 4945 




3.9 Mpc B 


3.95 


14.45 


359.3 


620.5 


10.2 


0.04 (C) 


HII/LI 


x,y,z,aa 


NGC 5253 




4.0 Mpc c 


2.61 


11.96 


30.51 


29.36 


9.3 


0.39 (W) 


HII (WR) 





Note. — Col.(l): Object name. ULIRGs arc listed first, and then LIRGs arc placed. The two very nearby galaxies, NGC 4945 and 5253 arc put last. Col. (2): 
Redshift. For the two very nearby sources, NGC 4945 and 5253, distances arc shown in Mpc. Col. (3)— (6): fi2, f25, f60 7 an d fioo are IRAS fluxes at 12 fim, 25 
//m, 60 /i.m, and 100 jixm, respectively, taken from Kim & Sanders (1998), Sanders ct al. (2003), or IRAS FSC catalog. Col. (7): Decimal logarithm of infrared 
(8—1000 fim) luminosity in units of solar luminosity (L Q ), calculated with Li R — 2.1 x 10 39 x D(Mpc) 2 x (13.48 x /i 2 + 5.16 x f 2 5 + 2.58 x / 60 + /ioo) ergs 
s _1 (Sanders & Mirabel 1996). Col. (8): IRAS 25-fim to 60-/J.m flux ratio. LIRGs with f25/feo < 0.2 and > 0.2 are classified as cool and warm sources (denoted 
as "C" and "W"), respectively, (Sanders ct al. 1988b). Col. (9): Optical spectral classification. "HII", "LI", "Sy2" , and "Unc" mean HH-rcgion, LINER, Seyfert 
2, and optically unclassified, respectively. For most sources, the classification is based on Vcillcux ct al. (1995) and Vcillcux ct al. (1999). For Superantennae 
and IRAS 20551 — 4250, the classification by Mirabel et al. (1991) and Due ct al. (1997), respectively, is adopted. For NGC 4945, we adopt the classification 
by Hcckman ct al. (1990) and Moorwood et al. (1996). For NGC 5253, the classification is from Conti (1991) and Walsh & Roy (1987), and "WR" means 
a Wolf Rayct galaxy. For NGC 7592, "W" and "E" indicate the western and eastern nuclei, respectively. Col. (10): Buried (obscured) AGN signatures and 
references, (a): Imanishi et al. (2007a). (b) Vcillcux ct al. (2009). (c): Nardini ct al. (2010). (d): Imanishi et al. (2006). (e): Tcng et al. (2008). (f): Dudley & 
Wynn- Williams (1997). (g): Imanishi & Dudley (2000). (h): Armus ct al. (2007). (i): Imanishi ct al. (2007b). (j): Imanishi ct al. (2008). (k): Imanishi et al. 
(2010b). (1): Risaliti et al. (2003). (m): Braito ct al. (2003). (n): Braito ct al. (2009). (o): Franccschini ct al. (2003). (p): Risaliti ct al. (2006). (q): Sani ct 
al. (2008). (r): Rousscl et al. (2006). (s): Imanishi (2006). (t): Imanishi ct al. (2009). (u): Maiolino ct al. (2003). (v): Evans ct al. (2008). (w): Lc Floc'h ct 
al. (2002). (x): Iwasawa et al. (1993). (y): Done ct al. (1996). (z): Guainazzi ct al. (2000). (aa): Itoh ct al. (2008). 

^Emission from the eastern (VV 114 E) and western (VV 114 W) nuclei combined. 
^Adopted from Bergman et al. (1992) and Maucrsbcrgcr ct al. (1997). 
C Adopted from Mauersberger ct al. (1997) and Thim ct al. (2003). 



Infrared 18 imaging of LIRGs 



TABLE 2 

Log for 18 /jm observations 



Object 


Date 


Integration time [min] 


Telescope/instrument 


(1) 


(2) 


(3) 


(4) 


Supcrantennac (1) 


2008 Sep 10 


15.5 


Gemini-S/T-ReCS 


IRAS 20551-4250 


2008 Sep 10 


15.5 


Gemini-S/T-ReCS 


IRAS 01004—2237 


2008 Sep 10 


31 


Gemmi-S/T-ReCS 


IRAS 04103-2838 


2008 Sep 10 


31 


Gemini-S/T-ReCS 


VV 114 E 


2008 Sep 10 


15.5 


Gemini-S/T-ReCS 


NGC 1377 


2008 Sep 10 


15.5 


Gemini-S/T-ReCS 


Supcrantennae (2) 


2009 Sep 8 


15.5 


Gemini-S/T-ReCS 


IRAS 00091-0738 


2009 Sep 8 


31 


Gemini-S/T-ReCS 


tt> A o nm oo riocr 

IKAb 00188— 0856 


2009 bep 8 


31 


Oemim-b/ 1-KeOb 


IRAS 01298-0744 


2009 Sep 8 


31 


Gemini-S/T-ReCS 


ESO 602-G025 


2009 Sep 8 


15.5 


Gemini-S/T-ReCS 


NGC 4945 


2010 May 31 


15.5 


Gemini-S/T-ReCS 


NGC 5253 


2010 May 31 


15.5 


Gemini-S/T-ReCS 


NGC 7592 


2010 Aug 30 


15.5 


Gemini-S/T-ReCS 


IRAS 08572+3915 (1) 


2008 Jan 24 


32 


Subaru/COMICS 


NGC 1614 


2008 Jan 24 


16 


Subaru/COMICS 


NGC 2623 


2008 Jan 24 


40 


Subaru/COMICS 


IRAS 08572+3915 (2) 


2009 Apr 7 


32 


Subaru/COMICS 


IRAS 15250+3609 


2009 Apr 7 


60 


Subaru/COMICS 


Arp 193 


2009 Apr 7 


56 


Subaru/COMICS 



Note. — Col.(l): Object name. The observed sources are sorted based on 
the telescope used and observing date. Sources above and below the horizontal 
solid line were observed with Gemini-South and Subaru, respectively. Superan- 
tennae and IRAS 08572+3915 were observed twice. Col. (2): Observing date in 
UT. Col. (3): Net on-source integration time in minutes. Col. (4): Telescope and 
instrument. 



TABLE 3 

Log for TV-band (8-13 fiu) observations, using Subaru COMICS 



Object 


Rcdshift 


Type 


Date 




Filter 


Integration 














Time [min] 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


I Zw 1 (PG 0050+124) 


0.061 


Syl a 


2006 Oct 


4 


Nil. 7 


26.7 


Mrk 1014 (PG 0157+001) 


0.164 


Syl a ' b 


2006 Oct 


4 


Nil. 7 


26.7 


PG 0844+349 


0.064 


Syl a 


2006 Oct 


4 


Nil. 7 


53.3 


PG 2130+099 


0.061 


Syl a 


2006 Oct 


4 


Nil. 7 


13.3 


PG 2214+139 


0.067 


Syl a 


2006 Oct 


4 


Nil. 7 


33.3 


Mrk 266 (NGC 5256) 


0.028 


Sy2 (N) + HII (S) b 


2009 Apr 


7 


N8.8 


13 



Note. — Col.(l): Object name. Col. (2): Redshift. Col. (3): Optical spectral classi- 
fication. a : Schmidt & Green (1983). b : Veilleux et al. (1999). For the LIRG Mrk 266, 
"N" and "S" denote the northern and southern nuclei, respectively. Col. (4): Observing 
date in UT. Col.(5): Filter used for the TV-band observation. TV8.8 (8.4-9.2 (im) or ATI. 7 
(11.2-12.1 /jm). Col. (6): Net on-source integration time in minutes. 
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TABLE 4 

Spitzer IRS OBSERVING LOG 



Object 


PID 


Date 


Intc 


:gration time 1 


sec] 






[UT] 


SL2 


SL1 


LL2 


LL1 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


IRAS 15250+3609 


105 


2004 Mar 4 


84 


84 


120 


120 


IRAS 20551-4250 


105 


2004 May 14 


56 


56 


56 


56 


NGC 1377 


159 


2005 Feb 18 


252 a 


252 a 


600 a 


600 a 


VV 114 E 


21 


2004 Jun 28 


210 b 


210 b 


180 b 


180 b 



Note. — Col.(l): Object name. Col.(2): PID number: 105 (PI = 
J. Houck), 159 (PI = R. Kennicutt), and 21 (PI = J. Houck). Col. (3): 
Observing date in UT. Col.(4)-(7): Net on-source integration time 
for SL2, SL1, LL2, and LL1 spectroscopy, respectively, in seconds. 



a IRS spectral mapping mode was employed. 

b IRS spectral mapping mode was employed. Total exposure time 
covering the E-nucleus was extracted. 



TABLE 5 

Flux measurement of LIRGs at ~18 im 



Object 


Flux (This work) 


Flux (Spitzer IRS) 


Standard star 


Std flux 




(mJy) 




(Jy) 


(1) 


(2) 


(3) 


(4) 


(5) 


Supcrantcnnae (1) 


595 


~600 


HR 7383 


2.4 


IRAS 20551-4250 


485 


~470 


HR 7869 


6.1 


IRAS 01004-2237 


330 


~270 


HR 400 


2.9 


IRAS 04103-2838 


290 


~250 


HR 1231 


34.6 


VV 114 E-l 


245 


~1000 (452+61 a ) 
~1000 (452+61 a ) 


HR 334 


6.1 


VV 114 E-2 


390 


HR 334 


6.1 


NGC 1377 


590 


~500 


HR 1231 


34.6 


Supcrantcnnae (2) 


640 


~600 


HR 7383 


2.4 


IRAS 00091-0738 


50 


~50 


HR 117 


2.6 


IRAS 00188-0856 


85 


~70 


HR 117 


2.6 


IRAS 01298-0744 


55 


~45 


HR 334 


6.1 


ESO 602-G025 


155 




HR 8015 


2.0 


NGC 5253 


5400 


~6000 (5500+539 b ) 


HR 5287 


7.3 


NGC 7592 W c 


210 




HR 8841 


2.6 


NGC 4945 A (1.8" x 1.8") 


200 


~1500 


HR 4906 


4.8 


NGC 4945 (3.6" x 8.9") 


1000 


~1500 


HR 4906 


4.8 


IRAS 08572+3915 (1) 


420 


~490 (552+46 d ) 


HR 3275 


8.7 


NGC 1614 


2970 


~2000 


HR 1231 


36.8 


NGC 2623 


260 


~300 


HR 3275 


8.7 


IRAS 08572+3915 (2) 


390 


~490 (552+46 d ) 


HR 3705 


25.1 


IRAS 15250+3609 


210 


~300 


HR 5763 


4.7 


Arp 193 


275 




HR 5219 


13.4 



Note. — Col.(l): Object name. The horizontal solid line separates sources observed 
with Gemini-South (above the line) and Subaru (below the line). Superantennae and IRAS 
08572+3915 were observed twice, and two flux measurements are shown separately. Col. (2): 
Flux measurement in this work in mJy. Col. (3): Spitzer IRS flux at 18 /im. References are 
Roussel et al. (2006), Brandl et al. (2006), Beirao et al. (2006), Imanishi et al. (2007a), Wu 
et al. (2009), and this paper (Figure 3). Col. (4): Standard star used for flux calibration. 
Col. (5): Adopted standard star flux at the observed wavelength in Jy. 

Measurement at 17.9 fim by Soifer et al. (2001) for VV 114 E. 

•^Measurement at 17.4 (5500+539 mJy; Moorwood & Glass 1982). A small 18.7 ^m 
flux measurement (2900 mJy) by Gorjian et al. (2001) was commented by Vanzi & Sauvage 
(2004). 

c Both the western (W) and eastern (E) nuclei were observed simultaneously, but the E 
nucleus was too faint to allow meaningful discussion. 
d Measurement by Soifer et al.(2000) at 17.7 /jm, using IRAS flux for a standard star. 



Infrared 18 imaging of LIRGs 



TABLE 6 

Emission size and surface brightness of LIRGs 



Ohiert 


PSF standard 




Physical size 


Surface brightness 


TL8 






l ell LOLL J 




CL^ kne^ 2 ! 




(1) 


(2) 


(3) 


(41 


(5) 


(6) 


Superantennae (1) 


HR 7383 


<0.14 


<160 


>1.0xl0 13 




IRAS 20551-4250 


HR 7869 


<0.15 


<120 


>0.72xl0 13 => >1.6xl0 13 


0.7 


IRAS 01004-2237 


HR 400 


<0.13 


<255 


>0.89xl0 13 




IRAS 04103-2838 


HR 1231 


<0.09 


<165 


>1.8xl0 13 




VV 114 E-l 


HR 334 


<0.19 


<70 


>0.20xl0 13 




VV 114 E-2 


HR 334 


0.78 


295 


~0.02xl0 13 




NGC 1377 


HR 1231 


<0.14 


<16 


>0.92xl0 13 => >2.5xl0 13 


0.85 


Superantennae (2) 


HR 7383 


<0.19 


<210 


>0.67xl0 13 




IRAS 00091-0736 


HR 117 


<0.18 


<350 


>0.07xl0 13 >0.24xl0 13 


1.05 


IRAS 00188-0856 


HR 117 


<0.10 


<215 


>0.39xl0 13 >0.75xl0 13 


0.55 


IRAS 01298-0744 


HR 334 


<0.18 


<400 


>0.08xl0 13 >0.30xl0 13 


1.1 


ESO 602-G025 


HR 8015 


0.4 


185 


~0.03xl0 13 




NGC 5253 


HR 5287 


0.32 


6.3 


~1.5xl0 13 




NGC 7592 W 


HR 8841 


<0.19 


<85 


>0.20xl0 13 




IRAS 08572+3915 (1) 


HR 3275 


<0.12 a 


<125 


>l.lxl0 i3 >3.3xl0 13 


0.9 


NGC 1614 


HR 1231 


1.9 


580 


~0.02xl0 13 




NGC 2623 


HR 3275 


<0.34 


<115 


>0.07xl0 13 




IRAS 08572+3915 (2) 


HR 3705 


<0.18 a 


<190 


>0.45xl0 13 => >1.3xl0 13 


0.9 


IRAS 15250+3609 


HR 5763 


<0.05 


<50 


>3.4xl0 13 >7.0xl0 13 


0.6 


Arp 193 


HR 5219 


1.0 


420 


~0.009xl0 13 





Note. — Col.(l): Object name. The horizontal solid line separates sources observed with Gemini-South 
(above the line) and Subaru (below the line). Col. (2): PSF standard star. Col. (3): Spatial extent of LIRG 18 
fim emission in apparent size in arcscc. Col. (4): Spatial extent of LIRG 18 /im emission in physical size in pc. 
Col. (5): Emission surface brightness in Lq kpc — 2 estimated from the 18 fim flux (luminosity) and measured 
size for the 18 /mi emission regions. Values after correction for the 18 fim flux-attenuation are also shown for 
LIRGs that display clearly detectable 18 (im silicate dust absorption features. No estimate is attempted for 
NGC 4945 because of its complex morphology. Col. (6): Optical depth of the 18 fim silicate dust absorption 
feature for detected sources, estimated based on our method (Imanishi ct al. 2007a; §4.2 of this paper). 

a Soifcr ct al. (2000) constrained to be <0'/22. 



TABLE 7 

Flux, size, and emission surface brightness of sources observed at ~10 /xm 



Object 


Flux (This work) 


Standard 


Std flux 


Size 


Physical size 


Surface brightness 




(mjy) 


star 


(Jy) 


(arcscc) 


(PC) 


(L Q kpc- 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


I Zwl 


380 


HR 106 


4.35 


<0.21 


<220 


>0.52xl0 13 => >2.6xl0 13 


Mrk 1014 


80 


HR 500 


5.15 


<0.14 


<360 


>0.33xl0 13 => >1.6xl0 13 


PG 0844+349 


50 


HR 3275 


20.3 


<0.10 


<110 


>0.29xl0 13 => >1.4xl0 13 


PG 2130+099 


135 


HR 8287 


3.90 


<0.28 


<300 


>0.10xl0 13 => >0.5xl0 13 


PG 2214+139 


75 


HR 8287 


3.90 


<0.17 


<190 


>0.16xl0 13 => >0.8xl0 13 


Mrk 266 N 


92 


HR 5219 


46.54 


0.56 


300 


0.02X10 13 => O.lxlO 13 


Mrk 266 S 


20 


HR 5219 


46.54 









Note. — Col.(l): Object name. Col. (2): Flux measurement at ~10 fxm (Mil. 7 or AT8.8) in mjy. Col. (3): Standard 
star used for flux and PSF calibration. Col. (4): Adopted standard star flux at the observed wavelength in Jy. Col. (5): 
Spatial extent at ~10 f-im. in apparent size in arcscc. Col. (6): Spatial extent at ~10 fim in physical size in pc. Col. (7): 
Emission surface brightness in L© kpc -2 estimated from the observed 11.7 fim (or 8.8 fim) flux (luminosity) and the 
physical size of emission regions. Values after correction for the infrarcd-unprobed emission (— fraction of energetic 
radiation unabsorbed by dust) arc also added for sources classified optically as Seyferts (§5.1.4). Mrk 266 S is barely 
detected at 10 fim., but is too faint to estimate the spatial extent and emission surface brightness in a reliable manner. 



14 



Imanishi et al. 



Supjerantennae 
* 

E 

UN 


HR 7383 
• 

E 


.1 RAS ■ 20651 -4250 
N 

e> . 


HR 7869 

N 

E> 




J RAS". 01 .004^237. 


HR 400 

n 

N 


IRAS 04103-2838 
■ *: 

N 

E^ 


HR 1231 
% 

N 

E^ 




Superantennae (2) 



HR 334 


' NGC WS^S 

... ... 


If- 

N 





HR 7383 
• 



E 





HR 117 




HR 4906 



N 



> 



P 



Fig. 1. — Qa (18.3 fim) images of LIRGs and corresponding standard stars, taken with Gemini-South T-ReCS. The image size is 8" X 
8", except VV 114 E and NGC 4945 (two larger figures) whose image sizes are 8" X 12" and 12" X 12", respectively, to include spatially 
extended structures. The size of each image is proportional to the actual ficld-of-view. The north (N) and east (E) directions are indicated. 
Upper and lower display levels are set arbitrarily and adjusted to individual sources, to make interesting emission structures (e.g., diffraction 
ring pattern) visually clear. The 18 fim emission of NGC 4945 is spatially extended. To better visualize the morphology, the image is 
Gaussian-smoothed by 2 pixels X 2 pixels. 
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Fig. 2. — Q17.7 (17.7 fim) images of LIRGs and corresponding standard stars, obtained with Subaru COMICS. The image size is 8" X 
8" for data taken in 2008, while it is 5" X 5" for data taken in 2009 (smaller figures in the bottom), because only 50 rows were read out, 
due to high background emission from Earth's atmosphere. The size of each image is proportional to the actual field-of-view. The upper 
and lower display levels are varied, according to individual sources, to show interesting emission patterns. The 18 /jm emission of Arp 193 
is spatially extended. To reduce noise and better visualize the morphology, the image is Gaussian-smoothed by 2 pixels X 2 pixels. 
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Fig. 3.— Spitzer IRS low-resolution spectra of IRAS 15250+3609, IRAS 20551-4250, NGC 1377, and VV 114 E. The dotted line 
represents the power-law continuum to measure the optical depths of the 18 fim silicate dust absorption features, determined from data 
points at A rcst = 14.2 /mi and 24 fim (Imanishi ct al. 2007a; Imanishi 2009; Imanishi ct al. 2010a). 



